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OBJECTIVE 


To investigate modeling, control, 
and computational issues associ- 
ated with elastic manipulators 


SCOPE 


Revolute joints 

Actuators at joints only 

Shadow beam approach 

Small elastic motion, and limit of 
such an assumption 

Nonlinear model 

Control issues 

Illustrative example 
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Pinned - free link 

Reference frame located at the pin 
joint; describes rigid-body motion. 
Elastic motion is defined with 
respect to this frame 




Notes: x is the position of the point in the undeformed configuration 
The beam rotates in a horizontal plane 
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u(x,t) is obtained by integrating 
d C = [ i - (-|^-) 2 1 1/2 dT i 

where 

r|(x,t) = x + s(x,t) 

C(x,t) = x - u(x,t) 

On integration, 

C(x,t) = T1 - l/2j o T| [ (-|^-) 2 ] do 
or 

u(x,t) = - s(x,t) + l/2f [(-^-) 2 ] do 

J do 

s(x,t) : axial vibration term 

Integral : results in centrifugal 

stiffening term 

Neglect axial vibration 
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Elastic Displacement 

u(x,t) = - u(x,t) k, + v(x,t) b 2 

Position 

r = £(x,t) b, + v(x,t) b 2 

Velocity 

f = - b,+ |7b2+ fiXr 

9t dt 

Kinetic Energy 

K = 1/2J r . r dm 

Potential Energy 

V = 1/2 J L EI(x) (v") 2 dx 

Lagrangian 

L = K - V 

Notes : (’) corresponds to partial derivative with respect to time, ( ) 
corresponds to spatial derivative, p is the mass per unit length, and EI(x) is the 


flexural rigidity 
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L = 1/2 J q L p (V 2 + 0 2 V 2 + 2 x v 0 + X 2 e 2 ) dx 

- 1/2 J L EI(v") 2 dx 

- 1/2 0 2 J Q L p x J o x (v0 2 da dx 

+ 1/2 O' 1/2 / 0 X ^ ((v,)2)da]2p dx 

+ 1/2 0 p v J 0 X '^■'(( V, ) 2 ) da dx 

- 1/2 0 J Q L p v J o x (v0 2 da dx 

+ 1/2 0 2 J Q L p/ 4 [/^(vO^a] 2 dx 

Assumed Modes 

N, 

v(x,t) = £ 4>iOO a i(t) 

i= 1 

^ i( x ) : Admissible functions 
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Define 


m ij = J 0 P^iCx) <t>j(x) dx 
ky = J q L EI(x) <j)j (x) (|)j (x) dx 


Sjj(x) = J 0 X <t>i(x) <J>j'(x) dx 
Pij = J 0 P X Sij(x) dx 

Sijkl = Jo S ij( x ) s kl( x ) dx 
qijk = J 0 P <t>k(x) Sij(x) dx 
r i= J 0 P x <t>i(x) dx 


Then, 

Jq X ( v ') 2 dx = Sij(x) aj aj 

J Q L p x [ J 0 X (v') 2 do ] dx = pjj aj aj 

JcX X 7r ((V,)2) d ° ]2p dx = 4 s *jkl a i a j a k a l 

| 0 l p v [ J 0 X d7 ((v,)2) da i dx = 2 qijk a i a j a k 
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Substituting in the Lagrangian, 

L = 1/2 I 0 0 2 + 1/2 my aj aj + 0 rj aj 

- l/2[kjj + (py - m^G 2 ] a s aj 

# » *0 

+ 1/2 Sy^j aj aj a k aj + 0/8 Sjj^j aj aj a^ aj 

* • 

® Qijk *4 — qjjk aj aj a^ 

Example : 

Beam parameters 

Cross-section : 6 in x 3/8 in 
Length = 3.6576 m (12 ft) 

p = 4.015 kg/m 

El = 756.65 N . m 

Admissible functions : Normalized 

eigenfunctions of a pinned-free beam 

m ij =S ij ; k ij =0) i 25 ij 

where 8y is the Kronecker delta 

Notes : The summation convention, XZ m ij a i a j = m ij a i a j> etc., will be em- 
ployed for conciseness - i.e., repeated indices in an expression indicate summa- 
tion over appropriate range. 
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Natural frequencies and Centrifugal stiffening coefficients 




r= i i 

j= 2 

1 

2 

3 

15.82 

51.282 

106.983 

6.397 

1.861 

-0.366 

1.861 

17.905 

6.195 


Coriolis terms, qjj k 


Other coupling terms s^ 

Note: S;in = s iik i = s iilk = s iilk 


-0.36' 

6.19! 

35.99 


k = 1 

k = 2 

k = 3 

-0.152 

0.415 

0.077 

0.143 

-0.144 

0.347 

0.008 

0.169 

-0.143 

0.415 

-0.175 

0.883 

-0.144 

0.152 

-0.196 

0.169 

-0.117 

0.145 

0.077 

0.883 

-0.178 

0.347 

-0.196 

0.171 

-0.143 

0.145 

-0.152 


1 = 1 

0.669 

0.099 

0.001 


= 2 


= 3 


1 

2 

3 









w?mu : itblied , bf iher® W the H terms 9ok an d s ijkI are small. In addition, 
?u! ,,P lie . d b y ‘ h 5 cub,c an d quartic powers of modal amplitudes V 
they will be dropped from further development. amplitudes, t 
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Retaining terms only up to quadratic in 
modal amplitudes, 


L= 1/2 [I 0 - (Pij - m 



- 1/2 kjj a; aj 

The equation for rigid-body motion is 

7-[(I 0 - (Pij- m ij ) aj ) 9 ] = T 
dt 

And the elastic motion is described by 

my aj + [ky + (Pij - my) 0 2 ] aj = T <))i(0), 

i = 1, 2, Ni 

Measurements at x = 0 

0, = 0 + v'(0,t) 

0j = 0 + v (0,t) 
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Choices for the control model 

1. Ignore elastic effects completely 

Control model : I 0 0 = T 

6 = 0 ] 0 = 0 , 

2. Rigid-body model, with the sha- 
dow frame angle properly 
extracted 

Control model : I 0 0 = T 
0=0,- v'(0,t) 

0=0,- v'(0,t) 


3. A few elastic modes are 
included, and the modal coordi- 
nates are approximated 

Control model: 

■57 Wo - (Pij - m ij) a, 3j ) 0 ] = T 

i,j < N 2 , N 2 < N, 

Notes : N, is the number of modeled modes. N 2 is the number 
of modes used for controller design. N, = 3 for the following 
simulation results. 6 
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4. Appropriate number of sensors 
used to obtain accurate modal 
coordinates. 


N 



1 


Control synthesis 
Computed torque method 
Pointwise-optimal control method 
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Open-loop Maneuver 


Fig. 1 : Torque Profile for Open-Loop Maneu 


Flexibility included 
Flexibility neglected 


Fig. 2 : Position Response of the Beam for the Torqu 








Feedback Control 



Actual velocity 

A ... 

_ . Desired velocity 

/ \ 

\ 


Fig. 3 : Rigid Model, Velocity response 



Actual velocity 
Oesired velocity 


Fig. 4 : One Flexible Mode Included in the Model 





Rodions/sec 
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Feedback Control, contd. 



Fig. 5 : Three Flexible Modes Included in the Model 



Fig. 6 : Comparison of Open- and Closed-Loop 
Torques 
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Effects of centrifugal stiffening 

1. Provides a strong coupling between 
the rigid-body and elastic motion 

2. Increases the stiffness of the struc- 
ture 

3. Reduces the effective rigid-body 
inertia term. Can cause it to vanish 
if the elastic motion is large. May 
have to modify the model, or vary 
the torque profiles. 
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Torque profiles and their Fourier 

coefficients 

Coefficients of 


Sin ('ftC-Vt)i Co 5 (y\ oJ„i) 

1 


Of\-1 

C-D 

o 

° \y * 

cn o0<, 


X! 1 

(i-c-fT) 

2 - 1 - ) 
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Computational Issues for Control of 

multi-link flexible robot arm 

1. The dynamic model can be arrived 
at by modeling each link indepen- 
dently and imposing constraints at 

the joints 

2 The link geometry may not be sim- 
ple 

3. s ijk „ q ijk , may not be negligible, and 
the control model may include all 
the terms 

4. The choice of admissible functions 
for each of the links may be 
different 

5. Sampling rates - should not excite 
elastic motion 

6. Control input computation may pose 
formidable burden. 
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The above issues can be adequately 
addressed by selecting pointwise- 
optimal control law for control input 
computations, where, the inputs can be 
computed at least one time step ahead. 
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Conclusions 

1. A complete model for control of a 
flexible link is developed 

2. Modeling issues are examined 
within the context of an example 

3. Several control issues are investi- 
gated 

4. It is shown that centrifugal stiffening 
effect on rigid-body motion is 
significant 

5. There is a strong coupling between 
rigid-body and elastic motions; 
ignoring this coupling results in 
gross inaccuracies in response. 
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